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We studied experimentally the effect of turbulent thermal diffusion in a multi-fan turbulence 
generator which produces a nearly homogeneous and isotropic flow with a small mean velocity. Using 
Particle Image Velocimetry and Image Processing techniques we showed that in a turbulent flow with 
an imposed mean vertical temperature gradient (stably stratified flow) particles accumulate in the 
regions with the mean temperature minimum. These experiments detected the effect of turbulent 
thermal diffusion in a multi-fan turbulence generator for relatively high Reynolds numbers. The 
experimental results are in compliance with the results of the previous experimental studies of 
turbulent thermal diffusion in oscillating grids turbulence (Buchholz et al. 2004; Eidelman et al. 
2004). We demonstrated that turbulent thermal diffusion is an universal phenomenon. It occurs 
independently of the method of turbulence generation, and the qualitative behavior of particle spatial 
distribution in these very different turbulent fiows is similar. Competition between turbulent fluxes 
caused by turbulent thermal diffusion and turbulent diffusion determines the formation of particle 
inhomogeneities. 

PACS numbers: 



I. INTRODUCTION 

The main goal of this study is to describe the ex- 
perimental investigation of the effect of turbulent ther- 
mal diffusion in a multi-fan turbulence generator. Tur- 
bulent thermal diffusion is associated with the correla- 
tion between temperature and velocity fluctuations in a 
turbulent flow with an imposed mean temperature gra- 
dient and causes a relatively strong non-diffusive mean 
flux of particles in the direction of the mean heat flux. 
This effect results in the formation of large-scale inhomo- 
geneities in particle spatial distribution whereby particles 
accumulate in the vicinity of the minimum of the mean 
fluid temperature. Turbulent thermal diffusion was pre- 
dicted theoretically by Elperin et al. (1996; 1997) and 
detected experimentally by Buchholz et al. (2004) and 
Eidelman et al. (2004) in oscillating grids turbulence. 

The mechanism of the phenomenon of turbulent ther- 
mal diffusion for inertial solid particles is as follows. In- 
ertia causes particles inside the turbulent eddies to drift 
out to the boundary regions between the eddies (i.e., re- 
gions with low vorticity and maximum of fluid pressure) . 
Therefore, particles accumulate in regions with maxi- 
mum pressure of the turbulent fluid. Similarly, there is 
an outflow of particles from regions with minimum pres- 
sure of fluid. In a homogeneous and isotropic turbulence 
without large-scale external gradients of temperature, a 
drift from regions with increased or decreased concentra- 
tion of particles by a turbulent flow of fluid is equiproba- 
ble in all directions, as well as pressure and temperature 
of the surrounding fluid do not correlate with the tur- 
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bulent velocity field. Therefore, only turbulent diffusion 
determines the turbulent flux of particles. 

In a turbulent fluid flow with a mean temperature gra- 
dient, the mean heat flux is not zero, i.e., the fluctuations 
of temperature and the velocity of the fluid are corre- 
lated. Fluctuations of temperature cause fluctuations of 
fluid pressure. These fluctuations result in fluctuations 
of the number density of particles. Indeed, an increase of 
pressure of the surrounding fluid is accompanied by an 
accumulation of particles due to their inertia. Therefore, 
the direction of the mean flux of particles coincides with 
that of the heat flux, and the mean flux of particles is 
directed to the region with minimum mean temperature, 
and the particles accumulate in this region (Elperin et 
al. 1996). 

The mechanism of turbulent thermal diffusion is as- 
sociated with a nonzero divergence of a particle velocity 
field. The latter is caused either by particle inertia or 
inhomogeneity of fluid density in a non-isothermal low- 
Mach number turbulent fluid flow. Therefore, the effect 
of turbulent thermal diffusion can be also observed in the 
suspension of non- inertial particles (e.g., particles with a 
Stokes time of the order of 10~® — f 0~^ s in air flows) or 
for gaseous admixtures in non-isothermal low-Mach num- 
ber turbulent fluid flows (see Elperin et al., 1997). Note 
that when we refer to compressible non-isothermal fluid 
flow with low Mach numbers, it means that div (p v) sa 0, 
where v is the fluid velocity and p is the fluid density. The 
latter implies that divv w — (v ■ ^)p/p ^ 0. In partic- 
ular, in a non- isothermal fluid flow with a temperature 
gradient divv -(v • 'V)p/p « (v • V)T/T ^ 0, where 
T is the fluid temperature. 

Numerical simulations, laboratory experiments and 
observations in the atmospheric turbulence revealed for- 
mation of long-living inhomogeneities in spatial distri- 
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bution of small incrtial particles and droplets in turbu- 
lent fluid flows (see, e.g., Wang and Maxcy 1993; Ko- 
rolcv and Mazin 1993; Eaton and Fessler 1994; Fessler 
et al. 1994; Maxey et al. 1996; Aliseda et al. 2002; 
Shaw 2003). The origin of these inhomogeneities was in- 
tensively studied by Elperin et al. (1996; 1997; 2000a; 
2000b). It was pointed out that the effect of turbulent 
thermal diffusion is important for understanding differ- 
ent atmospheric phenomena (e.g., atmospheric aerosols, 
smog formation, etc). In particular, the existence of a 
correlation between the appearance of temperature in- 
versions and the aerosol layers (pollutants) in the vicin- 
ity of the temperature inversions is well known (see, e.g., 
Csanady 1980; Seinfeld 1986; Flagan and Seinfeld 1988). 
Turbulent thermal diffusion can cause the formation of 
large-scale aerosol layers in the vicinity of temperature in- 
versions in atmospheric turbulence (Elperin et al. 2000a; 
2000b). Observations of the vertical distributions of pol- 
lutants in the atmosphere show that maximum concen- 
trations can occur within temperature inversion layers 
(see, e. g., Csanady 1980; Seinfeld 1986; Jaenicke 1987). 
The characteristic parameters of the atmospheric turbu- 
lent boundary layer are: the maximum scale of turbulent 
flow is L ~ 10"^ — 10^ cm; the turbulent fluid velocity in 
the scale L is u '--^ 30 — 100 cm/s; the Reynolds num- 
ber is Re = uL/v ~ 10^ — 10^ (see, e. g., Csanady 
1980; Seinfeld 1986; Blackadar 1997), where v is the kine- 
matic viscosity. For instance, for particles with material 
density Pp ~ 1 — 2 g / cm^ and radius a = 20 /um the 
characteristic time of formation of inhomogeneities is of 
the order of 1 hour for the temperature gradient li^/100 
m and 2 hours for the temperature gradient lif/200 m. 
The effect of turbulent thermal diffusion might be also of 
relevance in different industrial non-isothermal turbulent 
flows (Elperin et al. 1998). 

Turbulent thermal diffusion is a new and fundamen- 
tal phenomenon. Therefore, this phenomenon should be 
studied for different types of turbulence and different ex- 
perimental set-ups. Previously phenomenon of turbu- 
lent thermal diffusion was studied experimentally only in 
oscillating grids turbulence whereby the Reynolds num- 
bers were not so high (see for details Buchholz et al. 
2004; Eidclman et al. 2004). In order to study the ef- 
fect of turbulent thermal diffusion at higher Reynolds 
numbers we constructed a multi-fan turbulence gener- 
ator. Similar apparatus was used in the past in tur- 
bulent combustion studies (Birouk et al. 1996) and in 
studies of turbulence-induced preferential concentration 
of solid particles in microgravity conditions (Fallon and 
Rogers 2002). The multi-fan turbulence generator allows 
us to produce nearly homogeneous isotropic turbulent 
fluid flow with a small mean velocity. Using Particle 
Image Velocimetry and Image Processing Techniques we 
determined velocities and spatial distribution of tracer 
particles in the isothermal and non-isothermal turbulent 
flows. Our experiments with non- isothermal turbulent 
flows were performed for a stably stratifled fluid flow with 
a vertical mean temperature gradient which was formed 




FIG. 1: The scheme of the test section. 



by a cooled bottom wall and a heated top wall of the 
chamber. We found that particles accumulate in the 
vicinity of the bottom wall of the chamber (where the 
mean fluid temperature is minimum) due to the effect of 
turbulent thermal diffusion. Sedimentation of particles 
in a gravity fleld in these experiments was very slow in 
comparison with accumulation of particles caused by tur- 
bulent thermal diffusion. In the present study we demon- 
strated that turbulent thermal diffusion occurs indepen- 
dently of the method of turbulence generation. 



II. EXPERIMENTAL SET-UP 

A multi-fan turbulence generator includes eight fans 
(120 mm in outer diameter and with controlled rotation 
frequency of up to 2800 rpm) mounted in the corners of 
a cubic Perspex box and facing the center of the box. 
The Perspex box is a cube 400 x 400 x 400 mm^, with 
eight 272 mm equilateral triangles mounted in its corners 
used as solid bases for the fans (see Fig. 1). Each fan was 
calibrated separately, and the input current and rotation 
speed were measured and logged. We also tested one 
fan operating above an electric heat source. We placed 
a thermocouple in the vicinity of the fan's motor, and 
repeated this test for different temperatures (290, 310, 
330 K). These experiments showed that the fan rotation 
speed does not depend on the temperature. 

At the top and bottom walls of the Perspex box we 
installed two heat exchangers with rectangular 3 x 3 x 15 
mm^ flns. The upper wall was heated up to 343 K, the 
bottom wall was cooled to 283 K. Therefore, compera- 
tively large vertical mean temperature gradient (~ 92 
K/m) was formed in the core of the flow. The tem- 
perature was measured with a high-frequency response 
thermocouple (0.005 inches in diameter) which was glued 
externally to a wire (2.2 mm of Cu in diameter coated 
with 0.85 mm Teflon). The accuracy of the temperature 
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measurements was of the order of 0.1 K. We found that 
the temperature measurements affected the flow in a very 
small area around the wire. Two additional fans were in- 
stalled at the bottom and top walls of the chamber in or- 
der to produce a large mean temperature gradient in the 
core of the flow. Our measurements showed that these 
two additional fans only weakly affected considerably the 
homogeneity and isotropy of the turbulent flow. 

Velocity fields were measured using Particle Image Ve- 
locimetry (PIV) technique. The flow was seeded with 
incense smoke and was illuminated by a Surelite LSI-10 
(Contimmm) Nd:YAG pulsed laser with a power of 170 
mJ/pulse. The light sheet optics includes spherical and 
cylindrical Galilei telescopes with tuneable divergence 
and adjustable focus length. We used a progressive-scan 
12 bit digital CCD camera (pixel size 6.7/im x6.7/xm 
each) with a dual-frame-technique for cross-correlation 
processing of captured images. A programmable Timing 
Unit (PC interface card) generated sequences of pulses 
to control the laser, camera and data acquisition rate. 
The data was processed using standard cross correlation 
techniques (DaVis 7.0 code, LaVision, Gottingen). 

An incense smoke with sub-micron particles as a tracer 
for the PIV measurements was produced by high temper- 
ature sublimation of solid incense particles. Analysis of 
smoke particles using a microscope (Nikon, Epiphot with 
an amplification of 560) and a PM-300 portable laser par- 
ticulate analyzer showed that these particles have an ap- 
proximately spherical shape and that their mean diame- 
ter is of the order of 0.7/xm. In order to prevent from any 
thermal effects caused by the incense smoke generator, 
we placed it far away from the test section behind a wall, 
so that the incense smoke was transported through a five 
meter long pipeline and was cooled before it entered the 
test section. The smoke was feeded into the test section 
at room temperature. We measured the smoke tempera- 
ture inside the pipeline at two locations: 0.5 m and 3.5 m 
from the generator. At the first location the smoke tem- 
perature was 318 K, while at the second location it was 
294 K. The number density of smoke particles inserted 
to the test section in the experiments was of the order of 
10* cm-3. 

We determined mean and r.m.s. velocities, two-point 
correlation functions and an integral scale of turbulence 
from the measured velocity fields. A series of 130 pairs 
of images acquired with a frequency of 4 Hz were stored 
for calculating the velocity maps and for ensemble and 
spatial averaging of turbulence characteristics. The cen- 
ter of the measurement region coincides with the center 
of the chamber. We measured the velocity in flow area 
of 92 X 92 mm^ with a spatial resolution of 1024 x 1024 
pixels. These regions were analyzed with interrogation 
windows of 32 x 32 pixels. A velocity vector was de- 
termined in every interrogation window, allowing us to 
construct a velocity map comprising 32 x 32 vectors. The 
mean and r.m.s. velocities for each point of the velocity 
map (1024 points) were determined by averaging over 
130 independent maps, and then over 1024 points. The 




FIG. 2: The scheme of the temperature measurements. The 

wire that holds the thermocouple is inserted from the bottom 
of the chamber (some parts of the test section are not shown). 

Table 1 



Flow parameters in a multi-fan turbulence generator 



Horizontal and vertical directions 


Y 


Z 


Reynolds number Ke = uL/i/ 


703 


875 


Integral length scale L (mm) 


14.85 


16.4 


r.m.s. velocity u (m/s) 


0.71 


0.8 


Turbulouci' integral time sca.k' r = Lj u. (ms) 


20.9 


20.5 


Rate of dissipation e — u'' jL {rr? / s^) 


24.1 


31.2 


Taylor microscale \ = \/1^vt (mm) 


2.17 


2.15 


Kolmogorov length scale rj = LRe^^^** (/xm) 


109 


102 


Re\ = u\/v m. the Taylor microscale 


103 


115 



two-point correlation functions of the velocity field were 

determined for each point of the central part of the ve- 
locity map (16 X 16 vectors) by averaging over 130 inde- 
pendent velocity maps, and then over 256 points. Our 
tests showed that 130 image pairs contain enough data 
to obtain reliable statistical estimates. An integral scale 
L of turbulence was determined from the two-point cor- 
relation functions of the velocity field. For this end we 
used exponential approximation of the correlation func- 
tion since the experimentally measured correlation func- 
tion did not reach zero values. 

Spatial resolution of of the velocity measurements was 
about 2.9 mm for a probed area 92 x 92 mm^ with in- 
terrogation window 32 x 32 pixels. The maximum tracer 
particle displacement in the experiment was of the order 
of 8 pixels, i.e., 1/4 of the interrogation window. The av- 
erage displacement of tracer particles was of the order of 
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2.5 pixels. Therefore, the average accuracy of the velocity 
measurements was of the order of 4% for the accuracy of 
the correlation peak detection in the interrogation win- 
dow which was of the order of 0.1 pixels (see, e.g., Adrian 
1991; Westerweel 1997; 2000). 

In order to create large mean temperature gradient, 
the top and bottom fans were run at different speeds than 
the peripheral fans (see below). This regime was found 
empirically. Clearly, this introduces a weak anisotropy 
in velocity fluctuations in the vertical direction which is 
less than 10%. Note that the anisotropy of turbulence 
is not essential for the validation of the phenomenon of 
turbulent thermal diffusion. 

Fluid flow parameters (at the rotation speed 1500 rpm 
for eight corner fans and 2300 rpm for top and bottom 
fans) at the horizontal Y and vertical Z directions are 
presented in Table 1. In the experiments the maximum 
mean flow velocity was of the order of 0.1 — 0.2 m/s while 
the r.m.s. velocity was of the order of 0.7 — 1.1 m/s. 
Thus, the measured r.m.s. velocity was much higher 
than the mean fluid velocity in the core of the flow. Fig- 
ure 3 shows the two-point correlation functions of the 
velocity field at the horizontal and vertical directions. In 
particular, in Fig. 3 we plotted the longitudinal velocity 
correlation coefficients, f{y) = {uy{0)uy{y)) / {Uy{0)) and 
f{z) = {uz{0)uz{z)) / (!i^(0)), where u are the fluid veloc- 
ity fluctuations. Figure 3 and Table 1 demonstrate that 
the multi-fan turbulence generator produced a weakly 
anisotropic turbulent fluid flow with a small mean veloc- 
ity. 

The energy spectrum of the fluid flow is shown in 
Fig. 4. The one-dimensional longitudinal energy spec- 
trum was determined by a standard procedure. In par- 
ticular, we determined the Fourier components Uy{ky) 
and Uz{kz) of the fluctuating velocity field, and then de- 
termined (|wj,(fcy)p} and {\uz{kz)\'^) , where (...) is the 
ensemble averaging over 130 independent velocity maps. 
In Fig. 4 we plotted horizontal and vertical energy com- 
ponent spectra. As can be seen from Fig. 4, the energy 
spectrum is different from —5/3 law. However, our study 
showed that the existence of the effect of turbulent ther- 
mal diffusion is independent of the slope of the energy 
spectrum. 

Spatial particle number density distributions were ob- 
tained using a single frame from the double frame cap- 
tured for PIV measurements. For this purpose the in- 
tensity of laser light Mie scattering by tracer particles 
was recorded and averaged over 130 single frames. We 
probed the central 9.2 x 9.2 cm region in the chamber 
by determining the mean intensity of scattered light in 
32 X 16 interrogation windows with a size of 32 x 64 pix- 
els. The vertical distribution of the intensity of the scat- 
tered light was determined in 16 vertical strips, which are 
composed of 32 interrogation windows. Variations of the 
obtained vertical distributions between these strips were 
very small. We used spatial average across the strips and 
ensemble average over 130 images of the vertical distri- 
butions of the intensity of scattered light. The turbu- 
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FIG. 3: The two-point longitudinal correlation functions of 
the velocity field at the horizontal direction (filled squares) 
and vertical direction (unfilled squares). 
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FIG. 4: The energy spectrum of the fluid flow. 



lent diffusion coeflicient in the test section was of 
the order of ~ 40 cm^ / s. The turbulent diffu- 
sion time T^^ = L'^/D^ ~ 3 seconds in the scale of core 
flow Ld = 10 cm. Thus the steady state particle spatial 
distribution is reached during several seconds, while the 
turbulence integral time scale r ~ 2 x 10~^ s. The mea- 
surements were started several minutes after the seed was 
inserted in the chamber. 



III. TURBULENT THERMAL DIFFUSION: 
THEORY AND EXPERIMENT 

Now let us discuss the effect of turbulent thermal dif- 
fusion and compare the theoretical predictions with the 
experimental results. The mean number density of par- 
ticles N advected by a turbulent fluid flow is given by 



dN 
'dt 



-div[7V(V + Veff)-D^ViV] = 0, (1) 
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FIG. 5: Mechanism of turbulent thermal diffusion of inertial 
particles. 



Veff = -T(UpdivUp) = -D^{1 + k) 



VT 



(2) 



where Z)^ — (t/3)(u^) is the turbulent diffusion coeffi- 
cient, T is the momentum relaxation time of the turbu- 
lent velocity field, u arc fluctuations of fluid velocity, V 
is the mean fluid velocity, Up are fluctuations of parti- 
cle velocity, T is the mean fluid temperature. Coefficient 
K depends on particle inertia (the particle size a), the 
Reynolds number and the mean fluid temperature. In 
Eq. Q we neglected a small molecular mean flux of par- 
ticles caused by molecular (Brownian) diffusion, molecu- 
lar thermal diffusion (or molecular thermophoresis) and 
small particle terminal fall velocity. Equation (^) was 
previously derived by different methods (see Elperin et 
al. 1996; 1997; 1998; 2000b; 2001; Pandya and Mashayek 
2002; Reeks 2005). 

For non-incrtial particles, their velocity coincides with 
fluid velocity v = V + u, and divv « — (v • 'V)p/ p k, 
(v • ''^)T/T, where p and T are the density and temper- 
ature of the fluid, and V = (v). Therefore, the effective 
velocity of non-inertial particles Voff = — T(udivu) is 
given by Veff = —D^ ( VT) /T. Here we used the equation 
of state for an ideal gas, and neglected small gradients 
of the mean fluid pressure. This effective velocity causes 
an additional turbulent flux of particles directed to the 
minimum of mean fluid temperature (phenomenon of tur- 
bulent thermal diffusion). Note that turbulent thermal 
diffusion for non-inertial particles is the purely kinematic 
effect. Indeed, the equation for the instantaneous mass 
concentration, C = uipn/p, of non-inertial particles in 



non-isothermal flow reads 

8C 1 

_ + (v V)C- -div(i^pVC) , (3) 

where irip is the particle mass and n is the instantaneous 
number density of particles. For very small molecular 
diffusion D, this equation reads 



dC 

dt 



-f (v- V)Cw0 



(4) 



which implies that the mass concentration C is conserved 
along the fluid particle trajectory. In homogeneous tur- 
bulence all trajectories are similar. Therefore, the num- 
ber density of non-inertial particles n oc p, i.e., the num- 
ber density of non-inertial particles behaves locally as the 
fluid density. In particular, the location of the maximum 
of the number density of non-inertial particles coincides 
with the location of the maximum of the fluid density, 
and vice versa. For very small mean fluid pressure gra- 
dients, Vp/p ~ — VT/T, where p is the mean fluid tem- 
perature. Therefore, the location of the maximum of the 
mean number density of non-inertial particles coincides 
with the location of the minimum of the mean fluid tem- 
perature, and vice versa. The equation for the mean 
number density of non-inertial particles TV is equivalent 
to the following equation for the mean mass concentra- 
tion C — m.p N / p oi non-inertial particles: 

fin _ _ 1 

^ + {y.V)C^-_Aiy[D,pVC). (5) 
ot p 

If one ignores the turbulent thermal diffusion term in 
Eq. ^ for the mean number density of non-inertial par- 
ticles TV, then the equation for the mean mass concentra- 
tion C of non-inertial particles has an incorrect form. 

For inertial particles, their velocity Vp depends on 
the velocity of the surrounding fluid v. In particu- 
lar, for a small Stokes time Vp « v — Tpdv/dt -I- 0(Tp) 
(see Maxey 1987), where Tp is the Stokes time. Us- 
ing the Navier-Stokes equation it can be shown that 
divvp w div V -f TpAP/p + 0{t^) (Elperin et al. 1996). 
The effective velocity (j^J for inertial particles is given by 
Voff = —Z3ya(VT')/T', where the coefficient a — 1+K{a), 
K{a) (X. Tp (X a? and a is the particle size. Therefore, 
the mean particle velocity is Vp = V -I- Voff. Turbulent 
thermal diffusion implies an additional non-diffusive tur- 
bulent flux of inertial particles to the minimum of mean 
fluid temperature (i.e., the additional turbulent flux of 
inertial particles in the direction of the turbulent heat 
flux). In order to demonstrate that the directions of the 
turbulent flux of inertial particles and the turbulent heat 
flux coincide, let us assume that the mean temperature 
T2 at point 2 is larger than the mean temperature Ti 
at point 1 (see Fig. 5). Consider two small control vol- 
umes " a" and " 6" located between these two points (see 
Fig. 5), and let the direction of the local turbulent veloc- 
ity at the control volume " a" at some instant be the same 
as the direction of the turbulent heat flux {vl6) (i.e., it is 
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FIG. 6: Vertical temperature profile. Here ^ is a dimension- 
less vertical coordinate measured in units of the height of the 
chamber, and Z — aX, the bottom of the chamber. 



directed to the point 1). Let the local turbulent veloc- 
ity at the control volume " &" be directed at this instant 
opposite to the turbulent heat flux (i.e., it is directed to 
the point 2). In a fluid flow with an imposed mean tem- 
perature gradient, pressure p and velocity u fluctuations 
are correlated, and regions with a higher level of pres- 
sure fluctuations have higher temperature and velocity 
fluctuations. Indeed, using equation of state of the ideal 
gas it can be easily shown that the fluctuations of the 
temperature 9 and pressure p at the control volumes " a" 
are positive, and at the control volume " 6" they are neg- 
ative. Therefore, the fluctuations of the particle number 
density n are positive in the control volume "a" (because 
inertial particles are locally accumulated in the vicinity 
of the maximum of pressure fluctuations), and they are 
negative at the control volume (because there is an 
outflow of inertial particles from regions with a low pres- 
sure). The mean flux of particles (un) is positive in the 
control volume "a" (i.e., it is directed to the point 1), and 
it is also positive at the control volume " b" (because both 
fluctuations of velocity and number density of particles 
are negative at the control volume "&"). Therefore, the 
mean flux of inertial particles (un) is directed, as is the 
turbulent heat flux {u9), towards the point 1. 

The contribution of turbulent thermal diffusion to the 
turbulent flux of particles is given by 



where D^k^ is the coefficient of turbulent thermal dif- 
fusion, — aN = (1 -I- k)N is the turbulent thermal 
diffusion ratio, and the coefficient a = k^/N is the spe- 
cific turbulent thermal diffusion ratio. 

Neglecting the term iVVoff in Eq. for the mean 
number density of particles, we arrive at simple diffu- 
sion equation: dN/dt = D^AN, where we neglected a 
small mean velocity V. The steady-state solution of this 
equation is A'^ = const, i.e., a uniform spatial distribu- 
tion of particles. On the other hand, our measurements 
in both, the multi-fan turbulence generator and an oscil- 
lating grids turbulence generator (Buchholz et al. 2004; 
Eidelman et al. 2004) demonstrate that the solution N — 
const is valid only for an isothermal turbulent flow. Let 
us take into account the effect of turbulent thermal diffu- 
sion in Eq. . Then the steady-state solution of Eq. 
reads: VN/N = -aVT/T, which yields 



N 
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t — = 

To 
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^-D^a^N = -D^k,^ , (6) 
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where Nq — N{T — Tq) and Tq is the reference mean 
temperature. 

Now let us discuss the measurements of mean tempera- 
ture and particle spatial distribution in the multi-fan tur- 
bulence generator, and compare the experimental results 
with the theoretical predictions. The mean temperature 
vertical proflle in the multi-fan turbulence generator is 
shown in Fig. 6. Here Z is a dimensionless vertical coor- 
dinate measured in the units of the height of the chamber, 
and Z = at the bottom of the chamber. 

Measurements performed using different concentra- 
tions of the incense smoke in the flow showed that the 
distribution of the scattered light intensity normalized 
by the average over the vertical coordinate light inten- 
sity, is independent of the mean particles number density 
in the isothermal flow. In order to characterize the spa- 
tial distribution of particle number density, N oc E'^ /E, 
in a non-isothermal flow, the distribution of the scat- 
tered light intensity E for the isothermal case was used 
to normalize the scattered light intensity E'^ obtained in 
a non-isothermal flow under the same conditions. The 
scattered light intensities E'^ and E in each experiment 
were normalized by corresponding scattered light inten- 
sities averaged over the vertical coordinate. The ratio 
E'^ /E of the normalized average distributions of the in- 
tensity of scattered light as a function of the normalized 
vertical coordinate Z is shown in Fig. 7. A typical nor- 
malized E^ /E image in YZ plane is shown in Fig. 8. 
Inspection of Figs. 7-8 demonstrates that particles are re- 
distributed in a turbulent flow with a mean temperature 
gradient, e.g., they accumulate in regions with minimum 
mean temperature (in the lower part of the chamber). 

In order to determine the specific turbulent thermal 
diffusion ratio, a, in Fig. 9 we plotted the normalized 
particle number density = N/Nq versus the normal- 
ized temperature difference Tj = (T — Tq)/To, where Tq is 
the reference temperature and iVo = N{T = Tq). Figure 9 
was plotted using the mean temperature vertical profile 
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FIG. 7: Ratio /E of the normalized average distributions * 
of the intensity of scattered hght versus the normalized verti- 
cal coordinate Z. FIG. 8: A typical normalized E'^ /E image in YZ plane. Here 

Y and Z are the horizontal and vertical coordinates. 



shown in Fig. 6. The normaUzed local mean tempera- 
tures [the relative temperature differences [T — Tq)/Tq\ 
in Fig. 9 correspond to the different locations inside the 
probed region. In particular, in Fig. 9 the location of the 
point with reference temperature Tq'is Z = Q (the lowest 
point of the probed region with a maximum N) . In these 
experiments we found that the coefficient a ~ 2.68. 

The specific turbulent thermal diffusion ratio a in the 
experiments with oscillating grids turbulence (Eidclman 
et al, 2004; Buchholz et al., 2004) was a = 1.29 - 1.87 
(depending on the frequency of the grid oscillations and 
on the direction of the imposed vertical mean temper- 
ature gradient), while in the multi-fan turbulence a = 
2.68. The latter value of the coefficient a is larger 
than that obtained in the experiments in oscillating grids 
turbulence, where the Reynolds numbers were smaller 
than those achieved in the multi-fan turbulence genera- 
tor. Therefore, we demonstrated that the specific turbu- 
lent thermal diffusion ratio a increases with increase of 
Reynolds number. Note also that the experiments with 
oscillating grids turbulence were performed with two di- 
rections of the imposed vertical mean temperature gra- 
dient (for stable and unstable stratifications). 

The specific turbulent thermal diffusion ratio, a = 
1 + K,{a), comprises two terms, the first one (which equals 
to 1) is independent of the particle size and the second 
term depends on the size of particles. In particular. 




5 10 15 T^xlO^ 

FIG. 9: Normalized particle number density Nz = N/No 
versus normalized temperature difference Tz = {T — To)/To. 

K,{a) (X Tp (X a^, where a is the particle size. For non- 
inertial particles, K{a) = 0, and a = 1. The deviation 
of the coefficient a in both experiments from a = 1 is 
caused by a small yet finite inertia of the particles and 
also by the dependence of coefficient k on the Reynolds 
numbers. The exact value of parameter a for inertial par- 
ticles cannot be found within the framework of the theory 
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of turbulent thermal diffusion (Elpcrin ct al. 1996; 1997; 
1998; 2000b; 2001) for the conditions of our experiments 
(i.e., for large mean temperature gradients). However, 
in the experiments performed for different ranges of pa- 
rameters and different directions of a mean temperature 
gradient, and in two different experimental set-ups, the 
coefficient a was more than 1, that agrees with the the- 
ory. Therefore, we demonstrated that turbulent thermal 
diffusion occurs independently of the method of turbu- 
lence generation. 

The size of the probed region did not affect our re- 
sults. The variability of the results obtained in different 
experiments was within 0.5%. This is caused by variabil- 
ity of optical conditions, light intensity variations in a 
light sheet and errors in light intensity detection. There- 
fore, it can be concluded that error in particle number 
density measurements is less than 0.5%. Note that the 
contribution of the mean flow to the spatial distribution 
of particles is negligibly small. In particular, the normal- 
ized distribution of the scattered light intensity measured 
in the different vertical strips in the regions where the 
mean flow velocity and the coefBcient of turbulent dif- 
fusion may vary, are practically identical (the difference 
being only about 1%). 

The effect of the gravitational settling of small particles 
(0.5 — 1 /Ltm) is negligibly small (the terminal fall velocity 
of these particles being less than 0.01 cm/s). Due to the 
effect of turbulent thermal diffusion, particles are redis- 
tributed in the vertical direction in the chamber: parti- 
cles accumulated in the lower part of the chamber, i.e., in 
regions with a minimum mean temperature. Some frac- 
tion of particles sticks to the fan propellers and chamber 
walls, and the total number of particles without feeding 
fresh smoke slowly decreases. The characteristic time of 
this decrease is about 15 minutes. However, the spatial 
distribution of the normalized number density of parti- 
cles does not change over time. It must be noted that 
the accuracy of the measurements in these experiments 
(^ 0.5%) is considerably higher than the magnitude of 
the observed effect (~ 5%). Therefore, our experiments 



detected the effect of turbulent thermal diffusion in the 
multi-fan turbulence generator. These results arc in com- 
pliance with the results of the previous experiments in 
oscillating grids turbulence (see Buchholz et al. 2004; 
Eidelman et al. 2004). 



IV. CONCLUSIONS 

We studied experimentally the effect of turbulent ther- 
mal diffusion in a multi-fan turbulence generator using 
Particle Image Velocimetry and Image Processing Tech- 
niques. In a turbulent flow with an imposed vertical 
mean temperature gradient (stably stratified flow) par- 
ticles accumulate in regions of minimum mean temper- 
ature. Therefore, our experiments detected the effect of 
turbulent thermal diffusion in the multi-fan turbulence 
generator, i.e., non-diffusive mean flux of particles in the 
direction of the mean heat flux. Turbulent thermal dif- 
fusion is an universal phenomenon. In particular, using 
two very different turbulent flows created by oscillating 
grids turbulence generator (Buchholz et al. 2004; Eidel- 
man et al. 2004) and multi-fan turbulence generator, we 
demonstrated that the qualitative behavior of particle 
spatial distribution in non-isothermal turbulent flow is 
similar. The same physics is responsible for formation of 
particle inhomogeneities, i.e., competition between tur- 
bulent fluxes caused by turbulent thermal diffusion and 
turbulent diffusion. 
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